Spin gaps and spin-flip energies in density-functional theory 
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Energy gaps are crucial aspects of the electronic structure of finite and extended systems. Whereas 
much is known about how to define and calculate charge gaps in density-functional theory (DFT), 
and about the relation between these gaps and derivative discontinuities of the exchange-correlation 
functional, much less is know about spin gaps. In this paper we give density-functional definitions 
of spin-conserving gaps, spin-flip gaps and the spin stiffness in terms of many-body energies and 
in terms of single-particle (Kohn-Sham) energies. Our definitions are as analogous as possible to 
those commonly made in the charge case, but important differences between spin and charge gaps 
emerge already on the single-particle level because unlike the fundamental charge gap spin gaps 
involve excited-state energies. Kohn-Sham and many-body spin gaps are predicted to differ, and the 
difference is related to derivative discontinuities that are similar to, but distinct from, those usually 
considered in the case of charge gaps. Both ensemble DFT and time-dependent DFT (TDDFT) 
can be used to calculate these spin discontinuities from a suitable functional. We illustrate our 
findings by evaluating our definitions for the Lithium atom, for which we calculate spin gaps and 
spin discontinuities by making use of near-exact Kohn-Sham eigenvalues and, independently, from 
the single-pole approximation to TDDFT. The many-body corrections to the Kohn-Sham spin gaps 
are found to be negative, i.e., single-particle calculations tend to overestimate spin gaps while they 
underestimate charge gaps. 

PACS numbers: 71.15.Mb, 31.15.ej, 31.15.ee, 72.25.Dc 



I. INTRODUCTION 

There is hardly any electronic property of a system 
that does not depend on whether there is an energy gap 
for charge excitations, or for particle addition and re- 
moval. Similarly, there is hardly any magnetic property 
of a system that does not depend in some way on whether 
there is an energy gap for flipping a spin, or for adding 
and removing spins from the system. 

The reliable calculation of charge gaps^ from first prin- 
ciples is nontrivial and still faces practical problems (rel- 
evant aspects are reviewed below), but at least con- 
ceptually it is clear how charge gaps are to be defined 
and quantified within modern electronic-structure meth- 
ods, such as density- functional theory (DFT)££& On the 
other hand, much less is is known about how to calculate, 
or even define, spin gaps. 

In the present paper we show how to define and cal- 
culate the spin gap in spin-DFT (SDFT), and predict 
that such calculations will encounter a spin-gap problem 
similar to the band-gap problem familiar from applica- 
tions of DFT to semiconductors or to strongly correlated 
systems. 

Section [HI of this paper is devoted to charge gaps. In 
Sec. Ill Al we recapitulate the conceptual difference be- 
tween fundamental gaps and excitation gaps. In Sec. Ill Bl 
we then recall the quantitative definition of the funda- 
mental gap and related quantities, such as the single- 
particle gap, and particle addition and removal energies. 
Section !!! Cl summarizes key aspects of the derivative dis- 
continuity, while Sec. Ill Dl describes the connection be- 



tween gaps and discontinuities within the framework of 
ensemble DFT. Although the final results of these sec- 
tions are well known, our treatment is different from the 
usual one in so far as we introduce many-body corrections 
to the gap and derivative discontinuities in completely 
independent ways, related only a posteriori via ensemble 
DFT. This way of proceeding is useful for performing the 
generalization to the spin case. 

For both the fundamental gap and the optical exci- 
tation gap, the gapped degree of freedom is related to 
particles: either particles are added to or removed from 
the system, or particles are excited to higher energy lev- 
els within the system under study. In ordinary atoms, 
molecules and solids, these particles are electrons, and 
the particle gaps of many-electron systems are a key 
property in determining the functionality of today's elec- 
tronic devices. 

The last decade has witnessed an enormous growth of 
interest in another type of system, and in devices result- 
ing from them, in which the key degree of freedom is the 
spin. In the resulting field of spintronics, and the develop- 
ment of spintronic devices, one is interested in controlling 
and manipulating the spin degrees of freedom indepen- 
dently of, or in addition to, the charge degrees of freedom. 
Here, the issue of the spin gap arises, and a number of 
questions for electronic-structure and many-body theory 
appear: What is the energy required to add a spin to 
the system? What is the energy cost of flipping a spin? 
How do these concepts differ from the fundamental and 
optical gaps involving particles? Can we calculate spin 
gaps from spin-density-functional theory, and if yes, what 
type of exchange-correlation (xc) functional is required? 
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In Sec. IIII1 we answer these questions. 

In Sec. IIII Al we contrast spin gaps with charge gaps, 
and in Sec. IIII Bl we propose a set of many-body and 
single-particle definitions for quantities related to the 
spin gap, such as spin-flip energies and the spin stiffness. 
We take care to ensure that all quantities appearing in 
our definitions can, in principle, be calculated from con- 
ventional SDFT or time-dependent SDFT (TDSDFT), 
and try to make the definitions in the spin case as analo- 
gous as possible to the charge case. However, this analogy 
can only be carried up to a certain point, and important 
differences between charge gaps and spin gaps emerge al- 
ready at this level. As a simple example, we consider, 
in Sec. IIII CI the Lithium atom, for which we confront 
calculated and experimental spin gaps. 

In Sec. IIII Dl we then use ensembles DFT to relate 
the spin gap to a derivative discontinuity that is simi- 
lar to, but distinct from, the one usually considered in 
the charge case. Finally, in Sec. IIII Fl we investigate the 
connection to excitation gaps calculated from TDSDFT. 
Equations are given that allow one to extract the var- 
ious spin gaps and related quantities from noncollinear 
TDSDFT calculations. For illustrative purposes we eval- 
uate these for the Lithium atom, and compare the gaps 
and discontinuities obtained from time-dependent DFT 
to those obtained in Sec. IIII CI from time-independent 
considerations. 

Sec. IIVI contains our conclusions. 



II. CHARGE GAP 
A. Fundamental gaps vs. excitation gaps 

To provide the background for this investigation, let us 
first briefly recapitulate pertinent aspects of charge (or 
particle 1 ) gaps. 

While by definition all gaps involve energy differences 
between a lower-lying state (in practice often the ground 
state) and a state of higher energy, important differences 
depend on how the extra energy is added to the system 
and what degrees of freedom absorb it. Therefore, differ- 
ent notions of gap are appropriate for different purposes. 
For processes in which particles are added to or removed 
from the system, which is subsequently allowed to relax 
to the ground state appropriate to the new particle num- 
ber, the key quantity is the fundamental gap (sometimes 
also called the quasiparticle gap) which is calculated from 
differences of ground-state energies of systems with dif- 
ferent particle number. As such, it is relevant for instance 
in transport phenomena and electron-transfer reactions. 
If energy is added by means of radiation, on the other 
hand, the particle number does not change, and the rele- 
vant gap is an excitation energy of the TV-particle system. 
This excitation gap (sometimes also called the optical 
gap), is relevant, e.g., in spectroscopy. 

In first-principles electronic-structure calculations, ex- 
citation gaps are today often calculated from time- 



dependent density- functional theory (TDDFT). Funda- 
mental gaps, on the other hand, involve ground-state en- 
ergies of systems with different particle numbers, and 
should thus, in principle, be accessible by means of static 
(ground-state) DFT. However, it is well known that com- 
mon approximations to DFT encounter difficulties in this 
regard. In semiconductors, for example, calculated fun- 
damental gaps are often greatly underestimated relative 
to experiment, and in strongly-correlated systems such as 
transition-metal oxides, gapped materials are frequently 
incorrectly predicted to be metallic, i.e., to have no gap 
at all. The resulting band-gap problem of DFT has been 
intensely studied for many decades. 

A major breakthrough in this field was the discov- 
ery of the derivative discontinuity of the exact exchange- 
correlation (xc) functional of DFT, which was shown to 
account for the difference between the gap obtained from 
solving the single-particle Kohn-Sham (KS) equations of 
DFT, and the true fundamental gap££i£ The problems 
occurring in practice for semiconductors and strongly- 
correlated systems are therefore attributed to the fact 
that common local and semilocal approximations to the 
exact xc functional do not have such a discontinuity. The 
development of DFT-based methods allowing to nonem- 
pirically predict the presence and size of gaps in many- 
electron systems continues to be a key issue of electronic- 
structure theory and computational materials science. 

B. Definition of fundamental charge gaps 

The fundamental charge gap E g is defined as the dif- 
ference 

E g (N) = I(N) - A(N), (1) 

where the electron affinity (energy gained by bringing in 
a particle from infinity) and ionization energy (energy it 
costs to remove a particle to infinity) are defined in terms 
of ground-state energies of the TV-particle system, as 

A(N) = E(N)-E(N + l) (2) 
/(TV) = E(N-1)-E(N). (3) 

The order of terms in these differences is the conventional 
choice. The definition of the fundamental gap is in terms 
of processes involving addition and removal of charge and 
spin. The change in the respective quantum numbers is 
±1 in TV, and ±1/2 in S. In chemistry^ the average of 
/ and A is identified with the electronegativity of the 
TV-particle system: (/(TV) + A(N))/2 = x(TV). 

The corresponding Kohn-Sham gap is defined analo- 
gously as 

E 9 ,ks(N) = Iks(N) - A KS (N), (4) 

where I K S (N) = E KS (N -1) - E KS {N) and A KS (N) = 
E K s{N) - E KS (N + 1). Since the KS total energy is 
simply the sum of the KS eigenvalues, Eks = Sfc=i e k, 
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this reduces to Iks(N) = -£n{N) and A KS {N) = 
—€n+i(N), from which one obtains the usual form 

E g ,Ks(N) = e N+1 (N) - e N (N), (5) 

where en(N) and €n+i(N) are the highest occupied and 
the lowest unoccupied state of the iV-particle system, re- 
spectively. 

The fundamental gap can also be written in terms of 
KS eigenvalues by means of the ionization-potential theo- 
rem (sometimes known as Koopmans' theorem of DFT), 
which states 

I(N) = -e N (N) (6) 
A(N) = I(N + l) = -e N+1 (N+l), (7) 

so that I(N) ee Iks(N), and 

E g (N) = e N+1 (N+l)-e N (N). (8) 

Note that in contrast with the KS gap ([5]) these eigen- 
values pertain to different systems. 

The relation between both gaps is established by 
rewriting the fundamental gap as 



(9) 



which defines A xc as the xc correction to the single- 
particle gap. By making use of the previous relations 
we can cast A xc as ^ 10 ' 11 

A xc - e N+1 (N + l)-e N+1 (N) = A KS (N)-A(N). (10) 

The important thing to notice in these expressions is 
that, due to protection by Koopmans' theorem, the ion- 
ization energy does not contribute to the xc correction 
A xc , so that the correction of the affinity and of the fun- 
damental gap are one and the same quantity. Also, note 
that all of these definitions can be made without any 
recourse to ensemble DFT and without any mention of 
derivative discontinuities. 



C. Nonuniqueness and derivative discontinuities 



where T s [n] is the noninteracting kinetic energy func- 
tional, F[n] — T[n] + U[n] is the internal energy func- 
tional, expressed in terms of the interacting kinetic en- 
ergy T[n] and the interaction energy U[n], fi is the chem- 
ical potential, v ex t(r) is the external potential and v s (r) 
the effective KS potential. 

Both the effective and the external potential are only 
defined up to a constant, which does not change the form 
of the eigenfunctions. Consider now a gapped open sys- 
tem, connected to a particle reservoir with fixed chemical 
potential initially in the gap, and gradually change the 
constant. As long as the change is sufficiently small, the 
chemical potential remains in the gap, the density n(r) 
does not change, and the derivatives on the left-hand side 
of Eqs. (fT2|) and (fl3|) change continuously 

However, once the change in the constant is large 
enough to affect the number of occupied levels, the sit- 
uation changes: As soon as a new level falls below the 
chemical potential, or emerges above it, the number of 
particles in the system changes discontinuously by an in- 
teger, and the chemical potential adjusts itself to the new 
total particle number. For later convenience we call the 
two values of /x on the left and the right of integer particle 
number /i_ and fi + , respectively. 

When the right-hand side of Eqs. H2]) and ^ 
changes discontinuously, the left-hand side must also 
change discontinuously. This means that the functional 
derivatives of F[n] and T s [n] change discontinuously for 
variations 8n(r) such that N passes through an integer, 
and are not defined precisely at the integer. We can also 
argue conversely that if the functional derivatives existed 
at all n(r) they would determine the potentials uniquely. 
Since the potentials are unique only up to a constant, the 
functional derivatives cannot exist for the density varia- 
tions Sn arising from changing the potential by a con- 
stant. In a gapped system, these are the 5n integrating 
to an integer. 

Either way, we see that the indeterminacy of the po- 
tentials with respect to a constant implies that the func- 
tionals F[n] and T s [n] display derivative discontinuities 
for certain directions in density space along which the 
total particle number changes by an integer. This is the 
famous integer discontinuity of DFT^&l 



The basic Euler equations of DFT is^^ 
SE\n] 



Sn(r) 



(11) 



Since E[n] — F[n] + j d 3 r v ext (r)n(r) and Exs[n\ 
T s [n] + J d 3 r v s (r)n(r), this implies 



and 



5F[n] 
Sn(r) 

ST s [n] 
5n(r) 



/i - Vext(r) 



= /i- v s (r), 



(12) 



(13) 



D. Connection of discontinuities and gaps: 
ensemble DFT 



Up to this point we have defined A xc as a many-body 
correction to the single-particle gap, and deduced the 
existence of derivative discontinuities from noting the 
nonuniqueness of the external potentials with respect to 
a constant. These two conceptually distinct phenomena 
are related by ensemble DFT for systems with fractional 
particle number, describing open systems in contact with 
a particle reservoir £iiLL For such systems ensemble DFT 
guarantees that the ground-state energy as a function of 
particle number, E(N), is a set of straight lines connect- 
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ing values at integer particle numbers. 

For straight lines, the derivative at any A can be ob- 
tained from the values at the endpoints: 



-A = E(N+1)-E(N)= — 



SE 



and 

-J = E{N)-E(N- 



; on 



N+SN 



N-8N 



H+ 



M- 



(5n(r) 



5£ 



N+SN 

(14) 



i5n(r) 



The many-body fundamental gap is thus the derivative 
discontinuity of the total energy across densities integrat- 
ing to an integer: 



E„ 



I- A 



SE[r 



Sn(r) 



SE[r 



Sn(r) 



(16) 



N~SN 



N+SN 

This energy functional is commonly written as E = 
T s + V + Eh + E xc , where the external potential energy 
V and the Hartree energy Eh are manifestly continuous 
functionals of the density. Hence, the energy gap reduces 
to the sum of the discontinuity of the noninteracting ki- 
netic energy T s and that of the xc energy E xc . 

The entire argument up to this point can be repeated 
for a noninteracting system in external potential v s . The 
energy of this system is Eks — T s + V s , of which only the 
first term can be discontinuous. Hence the fundamental 
gap of the KS system is given by the discontinuity of T s 



E, 



g.KS 



ST s [n] 



Sn(r) 



ST s [n] 



N+SN 



Sn(r) 



(17) 



N~SN 



Returning now to the many-body gap, written as the 
sum of the discontinuities of T s and E xc , we arrive at 



E n 



Eg,KS + 



SE xc [n] 



Sn(r) 

or, by means of Eq. ©, 
8E xc \n] 



5E T 



N+SN 



Sn(r) 



N-SN 



A X c — 



Sn(r) 



SE X 



N+SN 



5n(r) 



(18) 



(19) 



N—SN 



This identifies the xc correction to the single-particle 
gap as the derivative discontinuity arising from the 
nonuniqueness of the potentials with respect to an ad- 
ditive constant £&tL 

Importantly, this connection is not required to define 
the xc corrections and neither is its existence enough to 
conclude that these corrections are nonzero. Many-body 
corrections to the single-particle gap can be defined inde- 
pendently of any particular property of the density func- 
tional (or even without using any density-functional the- 
ory), and whether for a given system these corrections 
are nonzero or not depends on the electronic structure of 
that particular system, and does not follow from the for- 
mal possibility of a derivative discontinuity, because this 
discontinuity itself might be zero. Thus, the question of 
the existence and size of xc corrections to the charge gap 
must be asked for each system anew. As we will see in 
the next section, the same is true for the spin gap. 



III. SPIN GAP 

We have provided the above rather detailed summary 
of the definition of the fundamental charge gap and its 
connection to nonuniqueness and to derivative disconti- 
nuities to prepare the ground for the following discussion 
of the spin gap. In order to arrive at a consistent DFT 
definition of spin gaps, we follow the steps outlined in 
the charge case: (i) define appropriate gaps and their 
xc corrections, (ii) use the nonuniqueness of the SDFT 
potentials to show the existence of spin derivative dis- 
continuities, and (iii) identify a suitable spin ensemble to 
connect the two. 



A. Spin gap vs. charge gap 

To introduce a spin gap or a spin-flip energy (see be- 
low for precise definitions) we consider processes in which 
only the total spin of the system is changed, while the 
particle number remains the same. There cannot be any 
definition in terms of particle addition and removal en- 
ergies, since in these processes the charge changes, too, 
which is not what one wants the spin gap to describe. In 
other words, the change of quantum numbers related to 
a spin flip is ±1 for the spin and for charge. Note that 
this is an excitation energy, where the excitation takes 
place under the constraints of constant particle number 
and change of total spin by one unit. This is the key 
difference to the previous section, from which all other 
differences follow. 



B. Definition of spin gaps: spin-flip energies and 
spin stiffness 

First, we define the spin up- flip energy and the spin 
down-flip energy in terms of many-body energies as 

E si+ (A) = E(N, S + 1) - E(N, S) (20) 
E sf -(N) = E{N,S-l)-E(N,S). (21) 

Here E(N, S) is the lowest energy in the A-particle spin- 
s' subspace, where S is the eigenvalue of the z-component 
of the total spin, and we assumed that spin-up and spin- 
down are good quantum numbers. This implies, in par- 
ticular, that spin-orbit coupling is excluded from our 
analysis. (Of course these definitions only apply if the 
respective flips are actually possible; in other words, if 
S does not yet have the maximal or minimal value for a 
given A.) 

The differences E s f + (N) and E s f~(N) are similar to 
the concepts of affinity and ionization energy, Eqs. @ 
and ([3]). However, affinities and ionization energies are 
always defined with the smaller value (of N) as the first 
term in the differences, whereas spin-flip energies are con- 
ventionally defined as final state minus initial state, i.e. 
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FIG. 1: Left: Spin-resolved single-particle (KS) density-of- 
states of a spin-polarized insulator. Two spin-flip energies 
and two spin-conserving gaps can be defined. Right: The 
half- metallic ferromagnet is a special case in which the gap in 
one spin channel (say spin up) is zero. In this case, there is 
only one spin-conserving gap, equal to the sum of both spin- 
flip energies, E s Jg + E 3 gg = E^ KS , and the KS charge gap 
is zero, due to the presence of the gapless spin down channel. 
Figure courtesy of Daniel Vieira. 



both spin-flip energies measure an energy cost. There- 
fore, the down-flip is the spin counterpart to the ioniza- 
tion energy, while the up-flip is the spin counterpart to 
minus the electron affinity. 

A more important difference is that the spin-flip en- 
ergies involve excited-state energies E(N, S + 1) and 
E(N, 5 — 1) of the iV-particle system, instead of ground- 
state energies, and in this sense are more similar to the 
optical gap in the charge case than to the quantities used 
in evaluating the fundamental gap. Alternatively, these 
energies can also be considered ground-state energies for 
sectors of Hilbert space restricted to a given total S, but 
we will not make use of this alternative interpretation in 
the following. 

KS spin-flip energies are related analogously to single- 
particle eigenvalues, according to 



,sf+ _ 

KS ~ 

sf- _ 

KS ~ 



e KT) - e H(i) 
e l(l) - e H(T)> 



(22) 
(23) 



where all energies are calculated at the same N, S. Here 
1(a) means the lowest unoccupied spin a state, and H (a) 
means the highest occupied spin a state. Similarly, L(a) 
and h(a) denote the lowest occupied spin a state and 
h(a) the highest unoccupied spin a state, respectively. 
This notation is nonstandard, but helpful, and further 
illustrated in Fig. Q] 

In the same way, we can also define the spin conserving 
(sc) single-particle gaps in each spin channel, as 



^g,KS 



E 



•sc,l 
g,KS 



e '(T) 
e Kl) 



(24) 
(25) 



The spin-conserving gaps and the spin-flip energies are 
necessarily related by E^g 



h g,KS 



E S J S 



E 



KS ■ 



look at Fig. [T] clarifies these definitions. If the system is 
non spin polarized, both spin-flip energies and both spin- 
conserving gaps become equal to the ordinary KS charge 
gap, which in our present notation reads e; — ejj. 

In the same way as for charge gaps, we can now also 
consider the sum and the difference of the spin-flip ener- 
gies. The sum 



E sf ~ + E sf+ 



(26) 



of the energies it costs to flip a spin up and a spin down 
is formally analogous to the fundamental gap ([T]), but 
with the important difference that E s involves excited- 
state energies. The unusual sign (sum instead of differ- 
ence) arises simply because both spin-flip energies mea- 
sure costs, whereas the affinity featured in Eq. |T]) mea- 
sures an energy gain. 

The formal analogy to Eq. ([I]) suggests that the quan- 
tity defined in Eq. (|2"6"]) be called the fundamental spin 
gap. In practice, however, the name spin gap is more 
appropriately applied to the individual spin-flip energies. 
The physical interpretation of their sum, Eq. (|26p , is re- 
vealed by expressing it in terms of the many-body ener- 
gies by means of Eqs. ([20]) and (|2"Tj) : 



E s = E(N, S + I) + E(N, S — 1) — 2E(N, S). (27) 

This is of the form of a discretized second derivative 
d 2 E(N,S)/dS 2 , which identifies E s as the discretized 
spin stiffness [we anticipated this interpretation when at- 
taching a subscript s for stiffness to the sum in Eq. ([26]) ]. 
We note that half of I— A is known in quantum chemistry 
as chemical hardness, which conveys a very similar idea 
as stiffness. Generically, we refer to all three quantities 
E s f~ , E s f + and E s as spin gaps. 

The spin electronegativity can be defined as half of the 
difference of the spin-flip energies, x s — [E s f~ — E s ^ + )/2. 
This quantity has the following interpretation: if x s > 0, 
it costs less energy to flip a spin up than to flip a spin 
down, whereas if X s < the down flip is energetically 
cheaper. 

The KS spin stiffness is defined as the sum of KS spin- 
flip energies, 



E s ,ks = Ev 



KS 



E 



KS ' 



or, with Eqs. (22]) and ([23]), 



E 



s,KS 



(28) 



(29) 



A 



This is analogous to ([5]) , except that in spin flips nothing 
is removed to infinity or brought in from infinity. Thus, 
differently from the KS ionization energy and electron 
affinity, the spin-flip energies require two single-particle 
energies for their definition instead of one, and in contrast 
with the KS charge gap the KS spin stiffness requires four 
single-particle energies instead of two. 

This missing analogy is physically meaningful: Con- 
ventional gaps are defined in terms of particle addition 
and removal processes and are ground-state properties. 
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TABLE I: Kohn-Sham energy eigenvalues (in eV) for the 
Lithium atom. The Is f, 2s j and Is j levels are oc- 
cupied. KS: energy eigenvalues obtained by inversion from 
quasi-exact densities. XX denotes exact exchange," and KLI 
is the Krieger-Li-Iafrate approximation.— 





KS° 


KS 6 


XX a 


KLI-XX 


LSDA 


— £isT 


55.97 


58.64 


55.94 


56.64 


51.02 


— £2sT 


5.39 


5.39 


5.34 


5.34 


3.16 




3.54 




3.48 


3.50 


1.34 


—tisi 


64.41 


64.41 


67.18 


67.14 


50.81 


— £2s| 


8.16 


5.87 


8.25 


8.23 


2.09 



°Ref. QJ 
^Ref. QJ 



TABLE II: Single-particle spin-flip energies (|34[1 and (|35[1 and 
spin stiffness ((29}, their experimental (Exp) counterparts, 
(|20p . (|21[l and (|26[) . and the resulting xc corrections defined 
in (|30|) , (|3ip and ((32}, for the Lithium atom. In the columns 
labelled KS we employ KS eigenvalues obtained from near- 
exact densities, while in the columns labelled XX, KLI-XX 
and LSDA we use approximate eigenvalues obtained from 
standard SDFT calculations. The experimental values were 
obtained using spectroscopic data for the lowest quartet state 
4 P° from Ref. as well as accurate wave-function based the- 
ory from Ref. 1 161 . All values are in eV. 





KS a 


KS 6 


XX a 


KLI-XX 


LSDA 


Exp 


E sf+ 


60.87 


60.87 


63.70 


63.64 


49.47 


57.41 


E sf- 


-2.77 


-0.48 


-2.91 


-2.89 


1.07 





E s 


58.10 


60.39 


60.79 


60.75 


50.54 


57.41 



To define pure spin gaps (i.e., spin-flip energies and spin 
stiffness) in which the charge does not change, we cannot 
make use of particle addition and removal processes but 
have to use spin flip processes instead. However, spin 
flips are excitation energies, and we must specify both 
initial and final states to define them properly. 

We also note that the many-body spin stiffness has no 
simple expression in terms of eigenvalues which would be 
analogous to Eq. ([8]). Such an expression would require 
the spin counterpart to Koopmans' theorem I(N) = 
Ifcs(N), which is not available for spin-flip energies. 
Hence, in general both spin-flip energies E s * ~~ and E s * + 
may be individually different from their KS counterparts 



E 



sf- 

KS 



and E 



«/+. 
ks ■ 



<f- = E *f 



KS 

sf+ _ E sf+ 
— ^KS 



A sf ~ 

XC 

A sf+ . 



(30) 
(31) 



We can, moreover, establish a relation between the many- 
body spin stiffness and the KS spin stiffness by rewriting 
the former as 



E x — E. 



s,KS 



E sf+ 
^KS 



A! 



(32) 



which defines A s xc as the xc correction to the KS spin 
stiffness. The important thing to notice in Eq. ([32]) is 
that there is no reason to attribute A xc only to the up- 
flip energy. This is a key difference to the charge case, 
where I(N) = Iks(N) and the xc correction could thus 
be attributed only to the electron affinity. Rather, the 
spin-flip corrections are connected by 

A S XC = A S J- + A S J+. (33) 



C. Example: the Li atom 

To give an explicit example of the quantities intro- 
duced in the previous section, we now consider the Li 
atom. For this system, KS eigenvalues £h] = £2s|, 
£ ;t = e 2p|' e U = £ 2si and ehi = eisi have been obtained 



A sf + 
Asf- 

A s 



-3.46 
2.77 
-0.69 



-3.46 
0.48 
-2.98 



-6.29 
2.91 
-3.38 



-6.23 
2.89 
-3.34 



7.94 
-1.07 
6.87 
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"Ref. 

6 Ref. [H, takin 



g -£2 P T 



3.54 eV from Ref. US 



by numerical inversion of the KS equation starting from 
near-exact densities (see Table [I]) . 

The KS spin-flip energies are obtained as 



sf+ _ 
ks — 



<sf- 
KS 



£2sl - £2sT- 



(34) 
(35) 



They are given in Table |TT1 together with the spin stiff- 
ness E s , see Eq. ([26]) . Table |TT] also presents the cor- 
responding experimental many-body energy differences 
for the Li atom, which were obtained using spectroscopic 
data for the lowest quartet state 4 P° and accurate wave- 
function based theoryj 15 ' 16 Relativistic effects and other 
small corrections included in the experimental data are 
ignored since they are too small on the scale of energies 
we are interested in. 

Table [H] also gives the xc corrections to the single- 
particle spin flip energies, see Eqs. (f3"0")) and (f3"Tj) . and 
the xc correction to the spin stiffness, A s xc , see Eq. (|32p . 
As a consistency test we verified that relation (|33|) , which 
connects the xc corrections of the spin-flip energies to the 
xc corrections of the spin stiffness, is satisfied. 

We also carried out calculations using the exact- 
exchange (XX) eigenvalues of Ref. [ID in order to sepa- 
rately assess the size of exchange and correlation effects. 
The resulting value of A s x — —3.38 eV indicates a larger 
(more negative) correction than in the calculation includ- 
ing correlation. An approximate KLI-XX calculation 14 
yields very similar results, while the LSDA data are com- 
pletely different and do not even reproduce the correct 
sign. 

Three of the required "exact" KS single-particle eigen- 
values are also reported in Ref. [TU (we use the result of 
Ref. [13 for the missing value of 62pt)- The value of t2 S { is 



7 



quite different than the value reported in Ref. LL2J (—5.87 
eV versus —8.16 eV), and consequently we obtain a rather 
different value of A* c (-2.98 eV versus -0.69 eV). Nev- 
ertheless, both sets of data sustain our main conclusions 
in this section: 

(i) Simple LSDA calculations give rise to serious qual- 
itative errors. As can be seen from Table UTI one obtains 
spin- flip energies that are drastically to small (E s f + ) or 
have the wrong sign (E s f~). The resulting xc corrections 
also suffer from having the wrong sign. These shortcom- 
ings of the LSDA are hardly surprising in view of its 
well-established failure to describe the charge gap. 

(ii) Even the precise KS eigenvalues do not predict the 
exact spin flip energies and spin stiffness, i.e. the xc cor- 
rections introduced in Scc. lIII Bl on purely formal grounds 
are indeed nonzero. The absolute size of these corrections 
implies that a simple KS eigenvalue calculation of spin 
gaps can be seriously in error. 

(iii) Exchange-only calculations overestimate (in mod- 
ulus) the size of the gap corrections. This implies that 
there is substantial cancellation between the exchange 
and the correlation contribution to the full correction. 
This is the same trend known for charge gaps. 

(iv) The xc corrections to both the up-flip energy and 
the spin stiffness turn out to be negative; in other words 
the KS calculation overestimates these quantities. This 
is the opposite of what occurs in the case of the funda- 
mental charge gap, which is underestimated by the KS 
calculation. We note that hints of an underestimation 
of the experimental spin-flip energies by KS eigenvalue 
differences have also been observed for half-metallic fer- 
romagnets. In the case of CrC>2, for example, Ref. [I?] 
reports experimental spin-flip energies in the range 0.06 
to 0.25 eV and compiles SDFT predictions that range 
from 0.2 to 0.7 eV (and in one case even 1.7 eV). 



D. Nonuniqueness and derivative discontinuities in 
SDFT 

Above we pointed out that the effective and external 
potentials of DFT are determined by the ground-state 
density up to an additive constant. However, this state- 
ment only holds when one formulates DFT exclusively 
in terms of the charge density, as we have done in dis- 
cussing charge gaps. It does not hold when one works 
with spin densities, as in SDFT, or current densities, as 
in current-DFT (CDFT). 

In these cases the densities still determine the wave 
function, but they do not uniquely determine the corre- 
sponding potentials. A first example of this nonunique- 
ness problem of generalized DFTs was already encoun- 
tered in early work on SDFT, for the single-particle KS 
Hamiltonianfi^ Later, this observation was extended to 
the SDFT many-body Hamiltonian ) 19 i 20 and further ex- 
amples were obtained in CDFT— and DFT on lattices^ 

Nonuniqueness is a generic feature of generalized (mul- 
tidensity) DFTs, consequences of which are still under 



investigationj 23 ' 24 i 25 ' 26 i 27 ' 28 In particular, Refs. [T^ and 
[20I already point out that the nonuniqueness of the po- 
tentials of SDFT implies that the SDFT functionals can 
have additional derivative discontinuities, because, if the 
functional derivatives of F and T s in multi-density DFTs 
such as SDFT and CDFT existed for all densities, they 
would determine the corresponding potentials uniquely. 
Very recently, Gal and collaborators 2 ^ pointed out that 
one-sided derivatives may still exist, and explored con- 
sequences of this for the DFT description of chemical 
reactivity indices. 

Just as in the charge case, derivative discontinuities re- 
sult from the nonuniqueness of the spin-dependent poten- 
tials, while corrections to single-particle gaps result from 
the auxiliary nature of Kohn-Sham eigenvalues. In the 
charge case, both distinct phenomena could be connected 
by means of ensemble DFT for systems of fractional par- 
ticle number. The question then arises if a similar con- 
nection can also be established in the spin case. This 
requires an investigation of spin-ensembles. 

E. Spin ensembles 

Consequences of the nonuniqueness of the potentials 
of SDFT for the calculation of spin gaps were already 
hinted at in Refs. [l^ and [2(| where it was pointed out 
that there may be a spin-gap problem in SDFT similarly 
to the well known band-gap problem of DFT. 

To make these hints more precise, we first recall, from 
the above, that the quantity usually called the spin gap 
is actually what we here called the spin-flip gap, and is 
analogous to the ionization energy or the electron affin- 
ity in the charge case, not to the fundamental particle 
gap. The spin-dependent quantity that is most analo- 
gous to the fundamental particle gap is the discretized 
spin stiffness of Eqs. (|26[) and (|27p . However, regardless 
of whether one focuses on the spin-flip energies or on the 
spin stiffness, the spin situation is not completely anal- 
ogous to the charge situation because both the spin-flip 
gaps and the spin stiffness are defined in terms of excited 
states of an iV-particle Hamiltonian, while charge gaps 
are defined in terms of ground-state energies of Hamilto- 
nians with different particle numbers. 

To identify a suitable ensemble, we write the energy 
associated with a generic ensemble of two systems, A 
and B, as 

E w = (l-w)E A +wE B , (36) 

where < w < 1 is the ensemble weight. If A and B have 
different particle numbers, Na and Nb = Na ± 1, this 
becomes the usual fractional-particle number ensemble, 
which is unsuitable for our present investigation where 
the involved systems differ in the spin but not the charge 
quantum numbers. 

A spin-dependent ensemble was recently constructed 
by Yang and collaborator s 29 ' 30 in order to understand the 
static correlation error of common density functionals. In 
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this spin ensemble, A and B have different (possibly frac- 
tional) spin, but are degenerate in energy. The constancy 
condition, whose importance and utility is stressed in 

Refs.mia arises directly from the restriction of the en- 
semble to degenerate states. While useful for the pur- 
poses of analyzing the static correlation error, this spin 
ensemble is too restrictive for our purposes, as it excludes 
the excited states involved in the definition of spin-flip 
gaps and of the spin stiffness. 

Ensembles involving excited states have been employed 
in DFT in connection with the calculation of excitation 
energies j 31 ' 32 i 33 Here A and B differ in energy but stem 
from the same Hamiltonian, with fixed particle number. 
Excited-state ensemble theory leads to a simple expres- 
sion relating the first excitation energy to a KS eigenvalue 
difference 32 



5EZ 



E B -E A = e 



M+l 



dw 



(37) 



where Eb and Ea are the energies of first excited and 
the ground state of the many-body system, respectively, 
and e%j are the highest occupied and lowest unoc- 



Z M+1 



cupied KS eigenvalues, and E™ c is the ensemble xc func- 
tional. Equation (|3"T|) holds for ensemble weights in the 
range < w < 1/2. Levy showed 3 -^ that the last term 
in this equation is related to a derivative discontinuity 
according to 



dEZ 



dw 



5EZ 



8E% 



6n ( r ) n =n^> Mr) 

(38) 

for w — * 0. Here n w = (1 — w)ua + wns is the ensemble 
density, and the discontinuity arises because even in the 
w — ► limit the ensemble density does contain an admix- 
ture of the state B with energy Eb > Ea and thus decays 
differently from no &s r —> oo£^ Levy developed his ar- 
gument explicitly only for the spin-unpolarized case, but 
already pointed out in the original paper that the results 
carry over to spin-polarized situations. 

In our case, we take A to be the ground state and B to 
be the lowest-lying state differing from it by a spin flip. 
To be specific, let us assume that the spin is flipped up. 
In this case we obtain from Eqs. (|3T[) and ([38]) in the 
limit w — > 0, and using our present notation, 



E sf+ (N) = E(N,S+1)-E(N,S) 



-KT) _ t H(i) 



^(T) _t H(l) 



dw 



n-" T 



(39) 
(40) 



se: 



w=0r 



l h n ll 



5n;(r) 



l h n l\ 



8ni(r) 



(41) 



E 



^ KS (N) + (42) 
for w — > 0. Equation (|42[) . which is the ensemble ver- 
sion of our Eq. (|3ip . illustrates that KS spin- flip excita- 
tions, too, acquire a many-body correction arising from 
a derivative discontinuity. 

In the particular case in which the spin flip costs no 
energy in the many-body and in the KS system, the pre- 
ceding equation reduces to Af/ jr„ = 0, which is the con- 
stancy condition derived in Refs. I^ t 30l for spin ensembles 
of degenerate states. 

We note that the KS eigenvalues and the discontinu- 
ity in Eqs. (|4"0"1) to (|42p must be evaluated by taking 
the w — > limit of the w-dependent quantities, while 
the quantities in Eq. (f3"Tj) have no ensemble dependence. 
This complicates the evaluation of spin-flip energies and 
their discontinuities, as defined in Sec. IIIIB1 from en- 
semble DFT. Therefore, we turn to still another density- 
functional approach to excited states in order to evaluate 
these quantities: TDDFT. 



F. Connection to TDDFT 

TDDFT has established itself as the method of choice 
for calculating excitation energies in atomic and molec- 
ular systems, and is making rapid progress in nanoscale 
systems and solids as welL 35 ' 36 In this section we will 
make a connection between the preceding discussion and 
TDDFT, which will allow us to derive simple approxima- 
tions for the xc corrections to the single-particle spin-flip 
excitation energies and the spin stiffness. 

To calculate the spin-conserving and the spin-flip exci- 
tation energies, it is necessary to use a noncollinear spin- 
density response theory, even if the system under study 
has a ground state with collinear spins (i.e., spin- up and 
-down with respect to the z axis are good quantum num- 
bers). In this way the spin-up and spin-down density 
responses can become coupled, and the description of 
spin-flip excitations (for instance, due to a transverse 
magnetic perturbation) becomes possible. In TDDFT, 
the spin-conserving and the spin-flip excitation energies 
can be obtained from the following eigenvalue equations, 
which are a generalization of the widely used Casida 
equations 3 ^ for systems with noncollinear spin^ 



<t<t' i' a' 



^i' i^a' a^aa.^a' ol' ^a' a' i' a ~r" K- 



act ,crcr 
iotaot'i' aa' a' 



x,, 



K. 



act ,(T a 
"iaaa'i' era' a' 



v cr.d a' 



-wXi, 



(43) 
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EE{ 



K 



io>aa',t'aa'a 



-'a cr i a 



K, 



a. ct,a a 
ia>aa.'i' aa' a' 



Y 



v cr, a (7 



= uY ir 



(44) 



where we use the standard convention that i,i and a, a' 
are indices of occupied and unoccupied KS orbitals, re- 
spectively, and a'a',acr' are spin indices, and Lu a 'a'i'a = 
fia'er' — a ■ Choosing the KS orbitals to be real, without 
loss of generality, we have 



K 



taaa.' .v aa' a 



•A") 



d A r / dV Via (r)Vw(r) 



x !^ aa , (r, r', ,)fr ff (r')4 v (r')(45) 

Here, the subscript indices of the matrix elements K refer 
to the KS orbitals in the integrand, and the superscript 
spin indices refer to the Hartree-irc kernel 

0, CTCT ,(r,r» = 5 -0^ + fZ^Ar,r',co) , (46) 

where the frequency-dependent xc kernel is defined as the 
Fourier transform of the time-dependent xc kernel 



f xc , ,( r t r' t') 



(47) 



n(r,t)=n°(r) 



Here, n(r,t) and n°(r) are the time-dependent and the 
ground-state 2x2 spin-density matrix, which follow from 
the DFT formalism for noncollinear spinsj 38 i 39 ' 40 i 41 ' 42 

Eqs. (|4"5|) and (|44[) give, in principle, the exact spin- 
conserving and spin- flip excitation energies of the system, 
provided the exact KS orbitals and energy eigenvalues are 
known, as well as the exact functional form of f^, aa , . 
We will now consider a simplified solution known as the 
single-pole approximation! 43 ! 44 It is obtained from the 
full system of equations (|43p , ([44]) by making the Tamm- 
Dancoff approximation (i.e., neglecting the off-diagonals) 
and focusing only on the H(<r) — > l(<r') excitations. In 
other words, we need to solve the 4x4 problem 



Y, 



Ha.la' 



luY, 



Ha.lo 



^2 [Sa'oi'SvcxWla'Ha + K Haia'Mala' 
aa' 

(48) 

For ground states with collinear spins, the only nonva- 
nishing elements of the Hartree-xc kernel are 



fHxc fHxc 

J TT,TT' JU,W 



cHxc 

r TT,ll' 



/ 



£XC fX 

/T1.TJ.5 J.\ 



(49) 



(notice that there is no Hartree term in the spin-flip 
channel), and the spin-conserving and spin- flip excita- 
tion channels decouple into two separate 2x2 problems. 
For the spin-conserving case, we have 



det 



cj tt 



■ U) 

M 



- u sc + M 



(50) 



where we abbreviate M aa > :Cra > = K^^,^ Hala ,(u!) and 
w<t'<t = (jJi<t',h<x = £ia' - £Ha- From this, we get the two 
spin-conserving excitation energies as 



-11 



M- 



TT.TT 



M 



± 



+ t( w TT - W U + M TT,TT ~ M U,llY 



1/2 



(51) 



with the spin-conserving Kohn-Sham single-particle gaps 
Egxs = The two spin-flip excitations follow 

immediately as 



E sf+ 
E sf~ 



where E s Jg = w n and E S j[ s = w iT . 

This gives a simple approximation for the xc correction 
to the spin stiffness: 



M- 
M 



11, Tl 



IT.J.T ) 



(52) 
(53) 



M 



M 



IT4T 



(54) 



Explicit expressions for f^, aa , can be obtained from 
the local spin-density approximation (LSD A), and we 
list them here for completeness (see also Wang and 
Ziegler^): 



£XC 

TTJT 

JXC 

f 



XC 

TT,U 



■m,u 



dn 2 
dn 2 
dn 2 



+ 2(1 
-2(1 
-2( 




■0 

dnd( 



d 2 e h xc 


. a-0 2 d 2 4 c 


dnd( 


n d( 2 


d 2 e h xc 


(1 + C) 2 d 2 c h xc 


dndC, 


n 8( 2 


(1- 


~C 2 )d 2 e h xc 



dC 2 



(55) 



where /? 



Ti,TT = fft,ll and /Tf.lT = AT,U' and ii: is under " 
stood that all expressions are multiplied by 5(r — r') and 

evaluated at the local ground-state density and spin po- 
larization, n (r) = nof(r) +n ^(r) and Co(r) = [«0T( r ) - 
noj.(r)]/no(r). For the xc energy density of the spin- 
polarized homogeneous electron gas we take the standard 
interpolation formula 



3 (»,C) 



3 (n,0) 



(i + C) 4/3 + (i - 



4/3 



2 4 / 3 



x [eS c (n,l)- e ^(n,0)]. 



(56) 



The case of exact exchange (XX) in linear response can 
be treated exactly, though with considerable technical 
and numerical effort . 45 i 46 A simplified expression of the 
XX xc kernel was developed by Petersilka et o/.— , and we 
have generalized their expression for the linear response 
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TABLE III: Top part: Spin-conserving and spin-flip excita- 
tion energies, calculated with LSDA and KLI-XX using dif- 
ferences of KS eigenvalues and TDDFT in the single-pole ap- 
proximation (|51[) - (|53[) . Bottom part: TDDFT xc corrections 
to the KS spin-flip excitation energies, from Eqs. (|52|l . (I53[l . 
and to the KS spin gap, Eq. ((54}. All numbers are in eV. 



LSDA 
KS TDDFT 



KLI-XX 
KS TDDFT 



Exact 
KS" Exp 6 



E sct 


1.83 


2.00 


1.84 


2.01 


1.85 1.85 


E scl 


48.72 


48.89 


58.90 


59.31 


56.25 56.36 


E sf + 


49.47 


48.23 


63.64 


62.12 


60.87 57.41 


E af ~ 


1.07 


0.99 


-2.89 


-2.97 


-2.77 0.0 


A sf + 




■1.24 




■1.52 


-3.46 


A s/ " 

'-^■xc 




0.08 




■0.07 


+2.77 


A" 




•1.32 




■1.59 


-0.69 



"Evaluated from the KS eigenvalues of Ref. 1 1 21 

'Spectroscopic data from Ref.E3 {E 3C ^^) and Refs-EHi (E 3f +) 



of the spin-density matrix. We obtain (details will be 
published elsewhere): 



N 



/fL t T (r, r ) — 



^ fcT (r)^ T (rQ^ T (r)^ T (rQ _ _ 

2^ tz zjTznz^nzR 



i,k 



|r — r'|n|(r)n-f(r') 
and similarly for ff, ,Jr, r'), and 



/?I,u(r,r') 



N- t ,N l 

E 

i , h 



- rV n T ( r ) n l ( r ) n T ( r ') n l ( r 

= /fT,lT( r » r ')- (58) 

Here, iVj and N± are the number of occupied spin-up and 
spin-down orbitals. 

We have evaluated Eqs. (|51 |l -(l54 |) for the spin-con- 
serving and spin-flip excitation energies of the Lithium 
atom involving the H(<j) and l(<r') orbitals. The LSDA 
and KLI-XX orbital eigenvalues that are needed as input 
are given in Table Q] 

The associated excitation energies are shown in Table 
IIII1 where we compare KS excitations, i.e., differences of 
KS eigenvalues, with TDDFT excitations obtained using 
the single-pole approximation described above. All in all, 
the TDDFT excitation energies are not much improved 
compared to the KS orbital eigenvalue differences. The 
main reason is that the LSDA and KLI-XX KS energy 
eigenvalues are not particularly close to the exact KS 
energy eigenvalues, and furthermore that the single-pole 
approximation is too simplistic for this open-shell atom. 

However, we observe that the xc correction A* c to the 
spin stiffness E s , when directly calculated within LSDA 
or KLI-XX using the TDDFT formula ([53]). is reasonably 
close to the exact value, and has the correct sign. This 
tells us that, even though the KS spin gap itself may be 
not very good, the simple TDDFT expression (j54"|) gives 
a reasonable approximation for the xc correction to it. 



IV. CONCLUSION 

The calculation of spin gaps and related quantities is 
important for phenomena like spin-flip excitations in fi- 
nite systems^ the magnetic and transport properties 
of extended systems such as half-metallic ferromagnets 17 
and, quite generally, in the emerging field of spintronics 
and spin-dependent transport. 

Our aim in this paper was to show how to define and 
calculate spin gaps and related quantities from density- 
functional theory. The proper definition of spin gaps in 
SDFT is by no means obvious, and the straightforward 
extrapolation of concepts and properties from the charge 
case to the spin case is fraught with dangers. There- 
fore, we started our investigation by disentangling two 
aspects of the gap problem that in the charge case are 
usually treated together: the derivative discontinuity and 
the many-body correction to single-particle gaps. 

On this background, we then provided a set of DFT- 
bascd definitions of quantities that are related to spin 
gaps, such as spin-conserving gaps, spin-flip gaps and 
the spin stiffness, pointing out in each case where pos- 
sible analogies to the charge case exist, and when these 
analogies break down. In particular, spin- flips involve 
excitations, while particle addition and removal involves 
ground-state energies. As a consequence, single-particle 
spin-flip energies involve two eigenvalues (and not one) 
and single-particle spin gaps involve four (and not two). 
Moreover, each spin- flip energy may have its own xc cor- 
rection (there is no Koopmans' theorem for spin flips). 

An evaluation of our definitions for the Lithium atom, 
making use of highly precise Kohn-Sham eigenvalues and 
spectroscopic data, shows that the many-body correction 
to spin gaps can indeed be nonzero. In fact, unlike what 
is common in the charge case, this correction turns out 
to be negative, i.e. the single-particle calculation overes- 
timates the spin gap while it underestimates the charge 
gap. While this result for a single atom is consistent 
with available data on half-metallic ferromagnets^ sim- 
ilar calculations must be performed for other systems be- 
fore broad trends can be identified. 

Next, we connected the many-body corrections to 
the spin gap and related quantities to ensemble DFT 
and to TDDFT. The former connection makes use of 
a suitable excited-state spin ensemble (different from 
the degenerate-state spin ensemble recently proposed by 
Yang and collaborators^^) and depends on a crucial in- 
sight of Levy2i regarding excited-state derivative discon- 
tinuities. The latter connection employs a noncollincar 
version of the Casida equations^ which we evaluate, 
again for the Lithium atom, within the single-pole ap- 
proximation, in LSDA and for exact exchange. 

The development of approximate density functionals 
and computational methodologies that permit the reli- 
able calculation of spin gaps and related quantities, in- 
cluding their many-body (xc) corrections, remains a chal- 
lenge for the future. 
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